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We show that the nonpcrturbativc quantum transport equations, the 'Kadanoff- 
Baym equations', can be be understood as the ensemble average over stochastic 
equations of Langevin type. For this we couple a free scalar boson quantum field 
to an environmental heat bath with some given temperature T. The inherent 
presence of noise and dissipation related by the fluctuation-dissipation theorem 
guarantees that the modes or particles become thermally populated on average in 
the long-time limit. This interpretation leads to a more intuitive physical picture 
of the process of thermalization and of the interpretation of the Kadanoff-Baym 
equations. 



1 Motivation 

Non-equilibrium many body theory had been traditionally a major topic of 
research for describing various (quantum) transport phenomena in plasma 
physics, in condensed matter physics and nuclear physics. Over the last years a 
lot of interest for non-equilibrium quantum field theory has now emerged also in 
particle physics. A very powerful diagrammatic tool is given by the 'Schwinger- 
Keldysh' or 'closed time path' (CTP) technique by means of non-equilibrium 
Green's functions for describing a quantum system also beyond thermal equi- 
librium. The resulting causal and nonperturbative equations of motion (by 
various approximations), the so called Kadanoff-Baym (KB) equations, have 
to be considered as an ensemble average over the initial density matrix charac- 
terizing the preparation of the initial state of the system. If the system behaves 
dissipatively, as a consequence of the famous fluctuation-dissipation theorem, 
there must exist fluctuations. The Kadanoff-Baym equations have then to be 
understood as an ensemble average over all the possible fluctuations. This in- 
herent stochastic aspect of the KB equations is what we want to point out and 
thus provide, as we believe, some new physical insight into its merely complex 
structure El. In what follows below we want to point out its intimate connection 
to Langevin like processes. 

As an elementary reminder of a Langevin process let us first briefly review 
the description of classical Brownian motion. Consider a heavy 'Brownian' 
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particle with mass M placed in a thermal environment obeying an effective 
Langevin equation, i.e. 

t 

Mi + 2 J dt'T(t - t') x(t')= £(*) . (1) 

— oo 

Here £(t) has to be interpreted as a 'noisy' source driving the fluctuations of 
the Brownian particle. For many applications is completely specified by a 
Gaussian distribution with zero mean and the correlation kernel / : 

I(t-t'):={(Z(m')}) =2TT(t-t'). (2) 

((...)) denotes the average over all possible realizations of the stochastic vari- 
able £(t). In fact, the simple relation between the dissipation kernel T and 
the strength I of the random force £(t) just stated is a manifestitaion of the 
fluctuation-dissipation theorem and (in the long time limit) is in accordance 

with the equipartition condition ((p 2 )) / (2M) = T/2. 

For a further physical motivation let us return to quantum field theory and 
already point out some similarities. One of the major present topics in quan- 
tum field theory at finite temperature or near thermal equilibrium concerns the 
evolution and behavior of the long wavelength modes. These modes often lie en- 
tirely in the non-perturbative regime. Therefore solutions of the classical field 
equations in Minkowski space have been widely used in recent years to describe 
long-distance properties of quantum fields that require a non-perturbative anal- 
ysis. A justification of the classical treatment of the long-distance dynamics 
of bosonic quantum fields at high temperature is based on the observation 
that the average thermal amplitude of low-momentum modes is large. For the 
low-momentum modes \p\ <C T (and for a weakly coupled quantum field the- 
ory) their (Bose) occupation number tib approaches the classical equipartition 
limit. The classical field equations should provide a good approximation for the 
dynamics of such highly occupied modes. However, in a correct semi-classical 
treatment of the soft modes the hard, i.e. thermal modes cannot simply be 
neglected, but it should incorporate their influence in a consistent way. In a 
recent paper □ it was shown how to construct an effective semi-classical action 
for describing not only the classical behavior of the long wavelength modes be- 
low some appropriate cutoff k c , but taking into account also perturbatively the 
interaction among the soft and hard modes. By integrating out the 'influence' 
of the hard modes on the two-loop level (for standard 4 -theory) the emerg- 
ing semi-classical equations of motion for the soft fields can be derived from 
an effective action and become stochastic equations of motion of generalized 
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Langevin typen, which resemble in their structure the analogous expression to 
(|l|) . The hard modes act as an environmental heat bath. They also guarantee 
that the soft modes become, on average, thermally populated with the same 
temperature as the heat bath. For the semi-classical regime, where \p\ <C T, 
one finds for the ensemble average of the squared amplitude 

^((\m\ 2 ))~^T*±-n B {E p ). (3) 

Such kind of Langevin description for the non-perturbative evolution of (super- 
)soft modes (on a scale of W\ ~ g 2 T <C T) in non Abelian gauge theories has 
recently been put forward 0. The understanding of the behavior of the soft 
modes is crucial e.g. for the issue of anomalous baryon number violation due 
to the diffusion of topological Chern-Simons charge in hot electroweak theory 
(see egU and references listed therein). 

In analogy to the Langevin description stated above we want to sketch in 
the following (pedagogical) study the effect of the heat bath on the evolution 
of the system degrees of freedom by means of the 'closed time path Green's 
function' (CTPGF) technique. For this we discuss a free scalar field theory 
interacting with a heat bath. 

2 Stochastic Interpretation of KB equations 

We start with the CTP action for a scalar bosonic field cj> coupled to an envi- 
ronmental heat bath of temperature T: 

m 2 )4>-] (4) 

4>- + <f>- E-+ 4> + + 4r £— 4>~] . 

The system starts to evolve from some initial density matrix. The interaction 
among the system and the heat bath is stated by an interaction kernel involving 
a self energy operator £ resulting effectively from integrating out the heat 
bath degrees of freedom. Schematically this is sketched in fig. [I]. In (Eh the 
self energy contribution from the heat bath is parametrized in the Keldysh 
notation by the four self energy parts, which can be expressed by means of the 
standard contributions £< and E > . Clearly, this self energy operator is the 
only quantity which might drive the system towards equilibrium. If the heat 
bath is at equilibrium, then the Kubo-Martin-Schwinger relation holds: 

£>(fc) = e ko/T t<{k). (5) 
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Figure 1: Schematic picture for the interaction of the system with the heat bath. 



(Here and in the following A(k) := A(ko, k) denotes the 4-dim. Fourier trans- 
form of a (stationary and translational invariant) quantity A(xi, x 2 ) — > A(x = 
x\ — x 2 ) with respect to the 4-dim. relative coordinate x.) 

From §Q) it is now straightforward to obtain the equations of motion for 
the characteristic two-point functions For the retarded propagator one has 

(-□ -m 2 - E rot )D rot = 5, (6) 

where E ret := 0(ti — t 2 ) [S > — S < ]. Additional dynamical information comes 
from the equation of motion of the propagator D < 

(-□-m 2 )D < -I] rct L> < -I] < J D av = 0. (7) 

This is just the famous KB equation. (||) and (0) determine the complete and 
causal (non-equilibrium) evolution for the two-point functions. 

To get more physical insight into the (effective) action (||) and in the 
equations of motion we introduce the following real valued quantities: 

s(xi,x 2 ) := -sgn(*i - t 2 ) (Z > {x 1 ,x 2 ) - T: < (x 1 ,x 2 )) = s(x 2 ,xi) , (8) 

a(xi,x 2 ) := - (T, > (xi,x 2 ) - S < (xi,a;2)) = -a(x 2 ,xi) , (9) 

I(x 1 ,x 2 ) := -— (£ > (xi,x 2 ) + S < (a;i,a;2)) = I(x 2 ,x 1 ) . (10) 

Our notion for s and a serves as a reminder for the respective symmetry prop- 
erties. It basically represents the standard decomposition of the real and imag- 
inary part of the Fourier transform of the retarded self energy operator E rct . 
s yields a (dynamical) mass shift for the 4> modes caused by the interaction 
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with the modes of the heat bath, while a is responsible for the damping, i.e. 
dissipation of the <j) fields. The important thing to point out will be that / 
characterizes the fluctuations. 

We first note that the CTP action (Q) can be written as 

d 4 x^ [cj) + (-□ - to 2 ) - 4>~ (-□ - to 2 ) <f)-] (11) 

d 4 Xl d 4 x 2 ^ [-(<f> + - 0") (s + a) (0+ + </T) + i (0+ - </T) I W + - <n] ■ 

This expression is identical to the so called influence functional given by Feyn- 
man and Vernon. To the exponential factor e lS in the path integral the (s + a) 
term contributes a phase while the / term causes an exponential damping and 
thus signals nonunitary evolution. 

The two relevant equations of motion are stated as 

(-a - m 2 - s - a)D Ict = 6, (12) 
(-□ - to 2 - s - a)D< + (a + U)D &Y = . (13) 

We see that the last equation is the only one where I occurs. 

For the interpretation of s, a and I consider the long-time behavior of these 
equations. In this case we can assume that the system becomes translational 
invariant in time and space and the boundary terms are no longer important. 
For the spectral function one immediately finds 

It becomes obvious that s = 3?S rct contributes an (energy dependent) mass 
shift while a = QYt rct causes the damping of propagating modes, a is related 
to the commonly used damping rate T via 

f(*)=<ffl. (15) 
k Q 

For D K one finds in the long-time limit the relation 

D < (k) = £> rot (fc)E<L) av (A:) 

= D Ict (k)[-d{k) ~U(k)]D av (k) = -2in(k)A(k), (16) 

where, by employing KMS condition (^), 

E < (fc) 1 
n(k) = £>(fc)-S<(fc) = S HBik0) ' (17) 
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which indeed shows that the phase space occupation number in the long-time 
limit becomes a Bose distribution with the temperature of the heat bath. 

It is now very illuminating to explicitly write down the relation between 
a(k) and I(k) using the definitions (§) and ©: 

In the high temperature (classical) limit (fco one gets 

I(k) = —2ia(k) , (19) 
fc 

or, employing |[|), 

I(k) = 2TT(k). (20) 

Recalling the discussion of Brownian motion in the introduction this compares 
favorably well with (||). The physical meaning of / as a 'noise' correlator will 
become obvious. The relation (^|) thus represents the generalized fluctuation- 
dissipation relation from a microscopic point of view by the various definitions 
of I, a and f through the parts E < and S > of the self energy. 

To see now more closely the connection to stochastic equations we decom- 
pose the influence action S as given in (^) in its real and imaginary part and 
write for the corresponding generating functional a 

z[j + ,r] 

:= fv\d> + ,(b-} p\(h + , 6-] e *Si<l> + .<p-]+i3 + <t> + +*r4>- 
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2?K] e-^^ 1 ? (v[<t> + ,<!>-] P [4> + ,4>-] e ^ + ^-]+;^ + +^-*-+i« 

= I J-nc - i/ • e,j i = ((^b ,+ +e,r-e]}) (21) 

with N := fD£ e - ^ 7 £ . The action entering the definition of Z' is no longer 
S, but only the real part of the influence action ([n]). The generating functional 
Z[j + ,j~] can thus be interpreted as a new stochastic generating functional 
Z'[j + + — £] averaged over a random Gaussian (noise) field £ with the 
width function /, i.e. 

«0»:=i Jvt-Oe-* 1 - 1 *. (22) 



6 



From the last definition we find that the (ensemble) average over the noise field 
vanishes, i.e. ((£)) = , while the noise correlator is given by 

(23) 

From the stochastic functional Z' a Langevin equation for a classical <p 
field can now readily be derived. Noting that the fields (</> + }{ on the upper 
branch and (4>~)^ on the lower branch are equal (and denoted as 0j in the 
following), its equation of motion derived from Z' takes the form 

(-□ -m 2 -s)(j) ( -afe = (24) 

This, indeed, represents a standard Langevin equation. The spatial Fourier 
transform of the Langevin equation (pjj) then takes the form 



^(k,t) + (m 2 + k 2 -2r(fc,A< = 0))^+2 / dt'T(k,t - t')^(k,t') = £(fc,i). 



(25) 

The analogy between this Langevin equation ( p5[ ) and the one for a single 
classical oscillator is obvious. The important difference, however, is the fact 
that the corresponding relations (|l^) and (||) between the respective noise 
kernel I and friction kernel T only agree in the high temperature limit. 



One can further ask to what extend the classical equations of motion (24) 



together with (23) are an approximation for the full quantum problem given by 
the equation of motion ( |l3] ) for D K . Inverting (|24|) one finds for the correlation 
function in the long-time limit 

- i({ (4>-)t)) = -iD™\(tf))D av = -iEF* IB™ . (26) 

Note that (^6|) is indeed the relation (||) advocated in the introduction to 
hold in the (semi-) classical regime. This one has to compare with the full 
quantum correlation function D < of (|l6l). One thus has that (—a — if) is 
approximated by —il. Of course this is justified, if \a\ <^ I holds. Using the 
microscopic quantum version (^|) of the fluctuation-dissipation theorem this is 
equivalent to coth (J^*) 3> 1. Thus in the high temperature limit or - turning 
the argument around - for low frequency modes, i.e. for k <C T, the classical 
solution yields a good approximation to the full quantum case. To be more 



precise: In simulations one has to solve the classical Langevin equation (24) 
and calculate n-point functions by averaging over the random sources. 

One can also write down the equations of motion for the quantum two- 
point functions with external noised by introducing 'noisy' two-point propaga- 



tors. Averaging the equation of motion over the noise fields according to (22) 
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one indeed rederives the KB equation. This demonstrates that the KB equa- 
tion can be interpreted as an ensemble average over fluctuating fields which are 
subject to noise, the latter being correlated by the sum of self energies S < and 
S > , i.e. from a transport theoretical point of view the sum of production and 
annihilation rate. We want to note once more that the 'noisy' or fluctuating 
part denoted by I inherent to the structure of the KB equation (Q) guarantees 
that the modes or particles become correctly (thermally) populated, as can be 
realized by inspecting (|l6|). 

We close our discussion by noting that one can also pursue to derive a 
standard kinetic transport equation for the (semi-classical) phase-space distri- 
bution f{x,k,t) including fluctuations S. The derived kinetic transport pro- 
cess has the structure of the phenomenologically inspired Boltzmann-Langevin 
equation. Our approach carried out inB has to be considered as a clear deriva- 
tion from first principles. Indeed it shows (nearly) a one to one correspondence 
to the phenomenologically introduced scheme. However, also some severe in- 
terpretational difficulties in the interpretation of the fluctuating phase-space 
density remain. We refer the interested reader to our discussion in El. 

3 Some further conclusions 

In our discussions we have elucidated on the stochastic aspects inherent to 
the (non-) equilibrium quantum transport equations. We have isolated a term 
denoted by / which solely characterizes the (thermal and quantum) fluctuations 
inherent to the underlying transport process. By introducing a stochastic 
generating functional the emerging stochastic equations of motion can then 
be seen as generalized (quantum) Langevin processes. What is changed, if we 
replace our toy model of a free system coupled to an external heat bath by a 
self-coupled and thus nonlinear closed system? In an interacting field theory 
of a closed system the KB equations formally have exactly the same structure 
as in our toy model. The important difference, however, is that the self energy 
operator is now described fully (within the appropriate approximative scheme) 
by the system variables, i.e. it is expressed as a convolution of various two-point 
functions. Hence, an underlying simple stochastic process, as in our case an 
external stochastic Gaussian process, cannot really be extracted. However, we 
emphasize that the emerging structure of the KB equations is identical. The 
decomposition of the self energy operator into its three physical parts (mass 
shift s, damping a, and fluctuation kernel /) can immediately be taken over. 
Hence these three parts keep their clear physical meaning also for a nonlinear 
closed system. 

One can also nicely demonstrate how so called pinch singularities^ are regu- 
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lated within the non-perturbative context of the thermalization process. These 
singularities do (and have to) appear in the perturbative evaluation of higher 
order diagrams within the CTP description of non-equilibrium quantum field 
theory. They are simply connected to the standard divergence in elementary 
scattering theory. The occurrence of pinch singularities signals the occurrence 
of (onshell) damping or dissipation. This necessitates in the description of the 
evolution of the system by means of non-perturbative transport equations. 

As a further application (discussed in the talk) for the direct use oLsemi- 
classical Langevin equations we want to mention the recent work in a: By 
applying a microscopically motivated Langevin description of the linear sigma 
model, one can investigate the stochastic evolution of a so called disoriented 
chiral condensate (DCC) in a rapidly expanding system, expected to occur in 
ultrarelativistic heavy ion collisions. Within such an approach one finds that an 
experimentally feasible DCC, if it does exist in nature, has to be a rare event, 
but still occuring with some finite and nonvanishing probability. The statistical 
distribution of final emitted pion number out of domains shows a striking 
nonpoissonian and nontrivial behaviour. One should indeed interpret those 
particular rarely occuring events as semi-classical 'pion bursts' similar to the 
mystique Centauro candidates. A further analysis of this unusual distribution 
by means of the cumulant expansion shows that the reduced higher order 
factorial cumulants exhibit an abnormal, exponentially increasing tendency 
and thus serves as a new and powerful signature. The occurence of a rapid 
chiral phase transition (and thus DCCs) might then probably only be identified 
experimentally by inspecting higher order facorial cumulants 9 m (m > 3) for 
taken distributions of low momentum pions. 
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